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Summary 


Stratospheric infrared solar absorption spectra recorded in the 
1300-2000-cm -1 region have ben analyzed to derive simultaneous sunset 
profiles of NO and N0 2 . The spectral data were obtained at sunset with the 
0.02-cirr 1 resolution University of Denver interferometer system during a 
balloon flight from Alamogordo, New Mexico, near 33° N on October 10, 1979. 
The technique of nonlinear least squares spectral curve fitting has been 
used to analyze the stratospheric spectra in regions of absorption by NO and 
N0 2 lines. To correct for the rapid diurnal changes in NO and N0 2 
concen-trations near sunset, a grid of time-dependent altitude profiles has 
been calculated for both gases with a detailed photochemical model. From 
these values, normalized diurnal correction factors have been derived for 
each scan and have been incorporated into the onion-peeling inversion 
procedure. The retrieved NO and N0 2 profiles are compared to previously 
published remote ijifrared and in situ measurements. The C0 2 profile was 
also determined from the analysis of the NO spectral region and is 
presented. 
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1. Introduction 

Retrievals of the stratospheric gas concentration profiles of NO and 
N0 2 from spectra recorded in the solar occultation mode are complicated by 
rapid changes at twilight which cause the concentration to be a function of 
solar zenith angle as well as altitude. Photochemical model calculations 
have shown that large errors in inferred NO concentration in the lower 
stratosphere will result at sunrise and sunset if this effect is neglected 
in the inversion procedure [Murcray et al., 1978; Boughner et al., 1980]. 

The magnitude of the error increases with decreasing altitude and can exceed 
100 percent for NO near 25 km [Boughner et al., 1980]. For N0 2 , the changes 
in concentration during sunrise and sunset are of a smaller magnitude, but 
the corresponding errors still amount to several percent in the lower 
stratosphere [Kerr et al., 1977; Murcray et al., 1978; Larsen and Boughner, 
1981]. 

Although in situ techniques avoid this problem and have been used to 
measure stratospheric concentrations of NO and N0 2 [cf. Patel et al., 1974; 
Maier et al ., 1978; Mihelcic et al., 1978; Horvath et al., 1983], solar 
absorption measurements are important because the simultaneous profiles of a 
large number of related species (e.g. H 2 0, CH 4 , 0 3 , HN0 3 , N 2 0) can be ob- 
tained frequently from the same data set [cf. Goldman et al ., 1980; 

Louisnard et al., 1983] in addition to pressure and temperature [Toth, 1977; 
Park, 1982; Rinsland et al ., 1983a]. Shortly, satellite solar occultation 
measurements of NO and N0 2 (as well as other species) will be made on a 
global scale by the Halogen Occultation Experiment [Russell et al ., 1977] 
and by the Atmospheric Trace Molecule Spectroscopy (ATMOS) experiment 
[Morse, 1980]. It is important, therefore, to develop techniques for 



improving the accuracies of NO and N0 2 inversions from solar absorption 
data. 

Boughner et al . [1980] have described a procedure for the calculation 
of normalized factors which can be incorporated into an onion-peeling 
retrieval to correct for diurnal changes during a solar occultation event. 
These factors are derived from time-dependent altitude profiles calculated 
with a photochemical model and express the variation in concentration with 
solar angle at a fixed altitude in terms of the 90° (tangent point) 
concentration. By normalizing, the need to know the absolute level of the 
profile is removed and instead it is determined from the retrieval process. 
Blatherwick et al . [1980] used this approach to derive a preliminary NO 
sunset profile from measurements of the equivalent widths of three isolated 
NO lines in solar absorption spectra recorded with the University of Denver 
balloon-borne interferometer system. In these calculations, the diurnal 
correction factors were estimated from NO distributions obtained with a 
relatively simple photochemical model [Murcray et al ., 1978]. The starting 
profiles for the photochemical calculations were taken from values in the 
literature, and a standard atmosphere pressure-temperature profile was 
assumed in both the photochemical calculations and in the retrieval 
analysi s. 

Although the diurnal model and the equivalent width method employed by 
Blatherwick et al . [1980] represent a reasonable and simple technique to 
obtain a first approximation to the NO sunset profile, it is important to 
know if further refinements in the photochemical calculations and in the 
analysis of the spectroscopic data would yield a significant improvement in 
the accuracy of the results for NO. For this reason, we have undertaken a 
reanalysis of the same set of University of Denver stratospheric spectra. 
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In addition to NO, the profile of N0 2 has also been retrieved. The detailed 
photochemical model of Cal 1 i s et al . [1983] has been used to calculate 
time-varying altitude profiles for both gases. These profiles have been 
integrated along the refracted ray to determine normalized diurnal correction 
factors for each of the spectra. To match the experimental conditions as 
close as possible, the profiles of related trace species obtained from the 
analysis of the same data set and correlative measurement of pressure and 
temperature have been used as inputs to these photochemical calculations. 

The stratospheric spectra have been analyzed with the technique of nonlinear 
least squares spectral curve fitting [cf. Chang et al., 1978; Niple, 1980; 
Park, 1983]. Although the reliability of this technique has yet to be 
studied in detail, it should at least in principle be more accurate than the 
equivalent width method because of its more complete use of the available 
information. An iterative photochemical model calculation/spectral inversion 
procedure has been adopted so that the NO and N0 2 profiles obtained from the 
analysis are consistent with the retrieval results. The final sunset NO and 
N0 2 profiles obtained from the analysis are compared with the previous values 
derived from the same set of spectra [Blatherwick et al., 1980; Niple et al., 

1980] and with other published infrared and in situ measurements. 

In addition, the profile of C0 2 has been derived from the analysis of 
the NO spectral region. Since there is considerable evidence that the C0 2 
mixing ratio is about 325 ppmv in the mid-stratosphere [cf. Volz et al., 

1981] , the retrieval of the simultaneous C0 2 profile provides a check for 
systematic errors in the inversion procedure. 

Research at Systems and Applied Sciences Corporation was supported by 
NASA. Research at the University of Denver was supported by NASA and NSF. 

We thank R. A. Toth of JPL for sending us calculated line parameters for the 
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6.2 ym N0 2 bands prior to publication. P. L. Rinsland of NASA Langley 
Research Center developed some of the programs used in the analysis. 

2. Photochemical Calculations 

A detailed time-dependent photochemical model [Callis et al., 1983] 
which incorporates all chemical reactions involving the H0 X , Cl x , N0 X , 
and 0 X families known to be important was used to calculate time-varying 
altitude profiles of NO and N0 2 for specified input distributions of temper- 
ature, pressure, and the long-lived species such as O 3 , H 2 0, CH 4 , N 2 0, and 
HNO 3 . Photolysis rates include single scattering, calculated numerically 
for a spherical atmosphere. Reflection of direct solar radiation from the 
surface and clouds was approximated with a surface albedo of 30 percent. 

Correction factors derived from two sets of photochemical calculations 
were used in the analysis. The initial calculation (case A) took profiles 
of the long-lived species obtained from a steady-state version of the 
photochemical model and temperature and pressure values from the 1976 U.S. 
Standard Atmosphere. The analysis of the spectral data with these factors 
provided a preliminary set of NO and N0 2 profiles for the second set of 
calculations (case B). Values for N 2 0s consistent with these NO and N0 2 
profiles were inferred from an approximate iterative scheme that equated the 
increase in N 2 0 5 at night to its daytime loss by photolysis. The equations 
describing the variations of N 2 0 5 at night allowed for the conversion of N0 2 
into N 2 0 5 , via N0 3 , and into C10N0 2 . A Cl x mixing ratio which becomes 
asymptotic to 2 ppbv in the upper stratosphere was adopted. Initially, no 
N 2 05 was assumed to be present at sunset. This amount was increased and 
adjusted iteratively until the calculated N0 X (N0+N0 2 ) concentrations at 
sunset matched the case A values. 
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Profiles of H 2 0, CH^, 0 3 , HN0 3 , and N 2 0 obtained from the same data set 
[Goldman et al., 1980, Rinsland et al 1982a; Rinsland et al 1984a] and 
correlative radiosonde and global satellite measurements of temperature and 
pressure (M. Gelman, National Meteorological Center, private communication, 
1981) were also used in the revised (case B) calculations. The diurnal 
calculations were made for solar zenith angles corresponding to the date and 
latitude of the measurements and were run for 2-1/2 days, starting at local 
noon, which was sufficient for the short-lived species to adjust to the 
concentrations of the long-lived species. After this time, the diurnal 
profiles were found to be reproducible to « 10 percent. Time steps were 
adjusted to limit truncation errors and ranged from microseconds near 
sunrise and sunset to « 1 hour near noon and midnight. Model calculated 
values of NO and N0 2 for case B agreed within 20-30 percent with the initial 
profiles (case A) retrieved for these species. 

3. Diurnal Correction Factors 


We next discuss the definition of the path-averaged diurnal correction 
factors used in the analysis. The geometry of a solar occultation event is 
shown schematically in figure 1 for three atmospheric paths with solar 
zenith angles 0 O , 0!, and 0 2 . The nomenclature presented here follows that 
used by Goldman and Saunders [1979]. Ray 0 corresponds to the atmospheric 
path for a spectrum recorded prior to local sunset (0 o <9O°); rays 1 and 2 
trace the paths for two successive low Sun scans with tangent altitudes Z i 
and Z 2 respectively. In our example, the layer boundaries are defined by 
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Fig. 1. Solar ray paths for occultation measurements obtained from float 
altitude Z. for high Sun (ray 0) and low Sun (rays 1 and 2) cases. 

The tangent altitudes of the low Sun scans are Z-j and Z^t respectively. 
The upper boundary of the atmosphere is Zq. 
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the top of the atmosphere Z 0 , the balloon float altitude Z b , and the 
tangent altitudes of the two low Sun scans. For simplicity, only a single 
layer has been shown for the range of altitudes above the balloon so that 
ray i passes through layer 0 and i inner layers. 

The infinite resolution transmittance of the atmosphere at frequency v 
corresponding to absorption by a gas along path i across layer 0 and i 
inner layers has been approximated by 

i 

T (i) = exp { - l [ - k (j) U(i ,j) ]} ( 1 ) 

V V 

j=0 

where k v is the monochromatic absorption coefficient (cm 2 /molecule) at v 
in layer j and U(i,j) is the total absorber amount (molecules/cm 2 ) along 
path i across layer j. The values of k v are calculated from the posi- 
tions, intensities, air-broadened halfwidths, and lower state energies of 
lines of the absorbing gas, assuming a Voigt line shape and the mass- 
weighted effective pressure and effective temperature for the path through 
the layer. The column amount is obtained from the integration of density 
along the refracted ray, 

U(i,j) » / p(z,e) ds ( z , 9) (2) 

path i 
layer j 

where ds(z,0) is the path increment, z is the altitude, and 0 is the local 
solar zenith angle. 
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According to the analysis of Boughner et al . [1980], the quantity 
U(i,j) can be rewritten as 

U(i,j) = / N(z,0) p(z, 90°) ds(z,0) (3) 

path i 
layer j 

where p(z, 90°) is the sunset density distribution of the gas and 

N(z,e) = p(z,0)/p(z, 90°) (4) 

is the normalized variation of the concentration with zenith angle relative 
to the sunset value at altitude z. The N(z,0) factor is obtained from the 
time-dependent photochemical calculations. The advantage of using this 
normalization is that it is not very sensitive to reasonable variations in 
the photochemical model inputs [Larsen and Boughner, 1981]. 

The path-averaged diurnal correction factor for path i across layer j, 
R(i,j), is defined as the ratio of the time-dependent column amount along 
the path, U(i,j), to the column amount that would be present with the sunset 
distribution p(z,90°), 

R(i,j) = U(i,j)/U s (i,j), (5) 

where 

U s (i»j) = / p(z» 90°) ds(z, 0) 

path i 
layer j 


( 6 ) 
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If R(i»j)>l» the time-dependent column amount along the path is larger than 
the corresponding sunset column amount; if R(i,j)<l, it is smaller. The 
path-averaged sunset density corresponding to path i through layer j is 
given by 


P s (i,j) = U s (1,j)/S(i # j), (7) 

where 

S(i.j) = / ds(z,e) (8) 

path i 
layer j 

is the total path length of the ray across the layer. 

In an onion-peeling inversion, an average density or average mixing 
ratio is inferred from the retrieved column amount in the tangent layer and 
then is usually assigned to this same altitude region in the analysis of 
spectra with successively lower tangent heights. However, if the 
concentration varies rapidly with altitude, different paths across a layer 
will correspond to different average densities (or mixing ratios), and this 
approximation will cause errors in the calculation of the column amounts. 

To correct for gradients in the concentration within a layer, a second set 
of normalized factors has been included in the analysis, defined by 

N(i,j) = P s (i,j)/P s (j,j), 


(9) 
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where p s (j,j) is the path-averaged sunset density for the case that layer 
j is the tangent layer. The values of R(i»j) and N(i,j) must be adjusted 
iteratively so that the p(z,90°) distribution used in the calculations is 
consistent with the shape of the retrieved vertical profile. 

From equations (1), (5), (7), (8), and (9), the infinite resolution 
transmittance at frequency v corresponding to absorption by a gas along 
path i across layer 0 and i inner layers is 


T v (i) = exp 


i-1 

I k v (j) R(i>j) N(i,j) S(i,j) p $ (j,j) 
j=0 


+ 'k v (i) R(i ,i ) S(i,i) P $ (i,i) . 


( 10 ) 


The first term in brackets contains the contribution to the total absorption 
from layer 0 above the balloon and the (i-1) inner layers above the tangent 
layer; the second term is the contribution of the absorption from within the 
tangent layer. For comparison with the solar absorption spectra, the 
monochromatic transmittances must be convolved with the effective instrument 
line shape function and scaled to simulate the wavelength dependence of the 
instrumental response and the solar flux. These instrumental factors have 
been determined in the nonlinear least squares analysis of the spectral 
data. The value of p s (i,i), the path-averaged sunset density in the 
tangent layer, is the only unknown in equation (10) and is also determined 
in the least-squares analysis. 



12 


4. Data and Analysis 

The spectral data were acquired at an apodized FWHM resolution of 
0.02 cm -1 with a Michel son-type Fourier transform interferometer (maximum 
path difference 50 cm) during a balloon flight from Alamogordo, New Mexico 
(32.8°N, 106. 0°W), on October 10, 1979. Sample stratospheric spectra 
obtained in both the NO and N0 2 regions have been published in an atlas 
format with line positions and identifications of the observed atmospheric 
and solar features [Goldman et al ., 1983]. The atlas also contains the same 
information for additional spectral regions recorded during the 1978 and 
1981 University of Denver balloon flights. The identifications given in the 
present work are taken from this atlas. The signal -to-rms noise is about 
100 in the NO and N0 2 regions. 

The application of the technique of nonlinear least squares spectral 
curve fitting to the analysis of University of Denver balloon-borne strato- 
spheric spectra is discussed in several papers [Goldman et al., 1980, Niple, 
1980; Niple et al., 1980; Rinsland et al., 1982a; Rinsland et al., 1983a] 
and, therefore, a description is not repeated here. In the present study, 
the NO and N0 2 profiles have been derived from six scans recorded at the 
float altitude of 33.0 ± 0.5 km. The correlative National Meteorological 
Center pressure-temperature profile assumed in the photochemical model 
calculations has also been adopted in the ray -tracing calculations with the 
FASATM program [Gallery et al ., 1983]. The National Meteorological Center 
measurements indicate a tropopause height of 14.9 km. The solar zenith 
angle corresponding to the time at which the zero path difference of the 
interferogram is recorded has been shown to be the appropriate effective 
solar zenith angle for airmass calculations [Park, 1982; Kyle and 
Blatherwick, 1984] and has been used for each scan in the ray -tracing 
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calculations. We note that in the previous analyses of NO and N0 2 
[Blatherwick et al., 1980; Niple et al., 1980], the mean solar zenith angles 
were used in the airmass calculations. 

Solar CO lines occur in both the NO and N0 2 spectral regions. Their 
absorption has been simulated with the Minnaert formula [cf., Kilston, 

1975], using the method described by Rinsland et al . [1982a]. The molecular 
constants reported by Guelachvili et al . [1983] were used to calculate the 
positions of these lines. 

4.1 NO Retrieval 

The 1903.0-1907.0 cm -1 and 1911.95-1917.00 cm" 1 intervals (464 and 470 
data points, respectively) were selected for the analysis. Each region con- 
tains four NO lines (actually unresolved lambda doublets separated by 
0.011 cm -1 ), a number of prominent atmospheric absorption lines of C0 2 and 
H 2 0, and numerous lines of solar CO. The NO lines are strongest in the 
92.41° and 93.36° scans and are weaker at larger and smaller zenith angles; 
the C0 2 and H 2 0 lines strengthen with increasing zenith angle over the 
entire range of solar zenith angles observed from float altitude 
(79.2-94.9°). The 93.78° and 94.14° scans also show several weak N 2 0 
lines. No additional species were identified as absorbers within these 
spectral intervals. 

The sunset density of NO and the mixing ratios of C0 2 and H 2 0 were 
included as unknowns in the nonlinear least-squares fit to each of the 
stratospheric spectra. The N 2 0 profile derived previously from analysis of 
the same set of stratospheric spectra was assumed [Rinsland et al., 1982a]. 
Weak channel spectra with periods of 0.24, 0.48, and 1.90 cm -1 were observed 
and modeled with the expressions derived by Niple et al . [1980]. 
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The spectral line parameters were taken mostly from the 1982 Air Force 
Geophysics Laboratory (AFGL) compilations [Rothman et al 1983a, b]. These 
listings include the NO positions and intensities generated recently by 
Gill is and Goldman [1982] and the water vapor line parameters of Flaud et al . 
[1981]. Improved CO 2 positions and intensities determined from an analysis 
of Kitt Peak laboratory spectra [Rinsland et al 1983b, c; Rinsland et al., ' 
1984b] were assumed. The positions of the stronger H 2 0 lines were updated to 
agree with the precise values measured by Guelachvili [1983]. Initial 
fittings suggested that the intensity of the H 2 16 0 line near 1904.761 cm -1 
(J'=4, K a '=4, K c '=l J"=4, K a "=l, K c "=4) is overestimated in the 
1982 AFGL compilation. A revised value was determined by comparison of the 
relative intensity reported by Guelachvili [1983] for this line on his 
spectrum 1556 with the relative intensities of several additional stronger 
and weaker lines for which Toth and Farmer [1975] have reported absolute 
intensity values. By interpolation, an absolute intensity at room 

1 04 

temperature for the 1904.761 cm - line of -1.8 x 10 cm/molecule 

(64 percent of the 1982 AFGL intensity) has been estimated; with this value, 

improved least-squares fits were obtained near this line for all spectra. 

The identified N 2 0 lines are R-branch transitions from the Vj + v 2 * band; 
their positions were taken as reported by Amiot and Guelachvili [1976]. 

Table I presents the NO densities retrieved with three successive 
refinements in the analysis. For each set of results, the NO densities 
obtained from fitting the two spectral intervals agreed to *» 10 percent and 
have been averaged. The sunset densities obtained assuming a non-time- 
dependent (spherically symmetric) NO distribution with constant NO density in 
a layer (that is R(i,j) = N(i,j) =1, for all paths and layers) are listed to 
show the results which would be obtained without corrections for diurnal 
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changes and gradients in the vertical profile. The case A and case B 
results were calculated with the same set of unknowns in the least-squares 
spectral analysis and assuming the R(i,j) and N(i,j) factors listed in 
Table II. Because the shape of the NO profile above the balloon float 
altitude cannot be determined from the data and has, therefore, been fixed 
by initial assumption, the average values listed for the layer above the 
balloon are enclosed in parentheses to indicate their greater uncertainty. 

An additional iteration of the photochemical model results was not made 
because, as discussed below, the differences between the case A and case B 
diurnal correction factors are relatively small and the measured NO profile 
has a rather large uncertainty, particularly below 25 km. 

The R(i,j) factors for cases A and B differ by less than 0.03, despite 
some rather large differences in the values assumed in the two sets of 
photochemical calculations. For example, to match the preliminary results 
for NO (and NO 2 ), the total odd nitrogen concentrations in the lower 
stratosphere were increased by nearly an order of magnitude for the case B 
calculations. The small sensitivity of normalized diurnal correction 
factors to perturbations in a number of model input parameters has been 
noted previously [Larsen and Boughner, 1981]. 

Table III compares the calculated case B diurnal correction factors and 
the values estimated by Blatherwick et al . [1980]. For the 85.85° and 
91.30° scans, the case B results indicate that diurnal corrections are 
negligibly small, in agreement with Blatherwick et al . [1980]. For the 
lower Sun scans, the present study yields slightly larger values of R(i,j) 
for the inner layers above the tangent layer and smaller values for the 
inner-most (tangent) layer, where Blatherwick et al . [1980] assumed R(i,j) 


to be 1.0. 
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The N(i,j) factors calculated for case A were obtained with a NO 
density profile which decreased by a factor of 7 from 31 km to 17 km. The 
preliminary (case A) retrieval, however, indicated a much smaller gradient 
in the NO profile, which is reflected in the smaller N(i,j) values 
calculated in case B for the layers below the balloon float altitude. 

For the low Sun scans, 30 percent - 47 percent of the total NO column 
amount results from the layer above the balloon. Since NO is likely to be 
distributed non-uni formly at high altitudes, the N(i,j) values assumed for 
this layer are important in the retrieval. For cases A and B, they have 
been calculated with the profile from the steady-state photochemical model. 
Above 55 km, we assumed the shape of the reference profile of Frederick 
et. al . [1983]. The calculated values range from 1.0 to 1.21. 

Recent in situ measurements of NO [Horvath et al., 1983] show 
considerable variability between 30 and 50 km at mid-latitudes, probably due 
to the action of transport processes. Our calculated N(i,j) factors for the 
layer above the balloon are, therefore, subject to an unknown amount of 
error. Because a high Sun scan of apparent angle e and a low Sun scan of 
apparent angle (18O°-0) have identical ray geometries above the balloon 
float altitude, the uncertainty in the retrieval for a particular tangent 
height can be minimized by using an appropriate high Sun scan in the onion- 
peeling analysis. For the lower stratosphere where the NO concentrations 
are small (solar zenith angles of « 94° for this experiment), the 
contribution to the total airmass from the path above the balloon float 
altitude is most important. Fortunately, the high-Sun scan available from 
this flight (85.85°) minimizes the uncertainties at these altitudes. In 
principle, N(i,j) factors for the layer above the balloon can be determined 
from measurements of average mixing ratios from a series of high Sun scans 
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covering a range of solar zenith angles. However, the NO lines are weak 
features in high Sun scans and measurement errors may make this procedure 
insufficiently precise. 

Figure 2 shows the observed and least-squares best-fit spectra for the 
1903.0-1907.0 cm" 1 interval of the 85.85° scan. The observed and calculated 
spectra are plotted in the lower panel; at top are the residuals (observed- 
calculated), expressed as a percentage of the maximum observed amplitude in 
the interval. The standard deviation of the residuals is 1.32 percent. 
Arrows mark the locations of the four NO lines; their absorption is 
relatively weak. Most of the observed features in the high Sun scan are 
solar CO lines, which have broader widths than the atmospheric lines. 

Figure 3 shows the results in the same region for the 93.36° scan 
(tangent altitude 22.3 km). The standard deviation of the residuals is 
1.38 percent. The standard deviation of the residuals for the other scans 
range from 1.34 percent to 1.42 percent. For all the scans, the inclusion 
of parameters to model the sinusoidal modulation of the background caused by 
channel spectra was critical to obtain an accurate fit in the region of the 
NO lines. 

Table IV compares the NO profile determined from the present study and 
the values derived from the same spectra by Blatherwick et al . [1980] with 
photochemical effects included in the analysis. The volume mixing ratios 
obtained in this work are 50-70 percent larger than retrieved by Blatherwick 
et al . [1980]. These differences are larger than can be accounted for by 
the different treatments of the photochemical effects. As can be seen from 
the fit to the background in figure 3, channel spectra occur in the NO 
spectral region, and it is possible that some of the differences in the 
retrieved profile result from the difficulty in defining the background for 
measurements with the equivalent width technique. 
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Fig. 2 


. Comparison of the 1903. 0-1907. 0-cm region of a solar absorption 
spectrum (solid line) obtained at an astronomical zenith angle 
of 85.85° and a least-squares best fit to the data. The arrows 
mark the locations of four NO absorption lines. Residuals 
expressed as a percentage of the maximum observed amplitude 
are shown in the upper panel. 
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Fig. 3. Comparison of the 1903. 0-1907. 0-cm~^ region of a solar 

absorption spectrum (solid line) obtained at an astronomical 
zenith angle of 93.36° (tangent altitude = 22.3 km) and a 
least-squares best fit to the data. The results are shown 
in the same format as Fig. 2. 
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Most of the published NO measurements were obtained in situ under high 
Sun conditions. Figure 4 compares our profile (converted from density to 
volume mixing ratio) and published in situ measurements for the same 
latitude and season. The measurements of Patel et al . [1974] were obtained 
with a spin flip Raman laser at 28 km on October 19, 1973. We have plotted 
their measured high Sun mixing ratio, which was obtained about 5 hours after 
visible sunrise. The other in situ data were derived from measurements made 
with chemiluminescence detectors. Loewenstein et al . [1978] reported NO 
concentrations as a function of latitude. The rectangles in figure 4 
represent the range of their fall measurement near 33°N. Profiles derived 
from three flights in the fall of 1977 and 1978 by Ridley and Schiff [1981] 
are also plotted as are the measurements of M. McGhan and co-workers (as 
quoted by Ridley and Schiff, 1981). 

Our values are about a factor of 2 lower than the in situ measurements 
presented in figure 4. This difference is expected since NO is rapidly 
converted to NOa near sunset. For example, the photochemical model 
calculations predict that the ratio of the sunset concentration to the 
concentration in the afternoon at a solar zenith angle of 45° is 0.56 at 
17.5 km, 0.61 at 23.5 km, and 0.73 at 31.0 km. Clearly, quantitative 
comparisons between results obtained with absorption and in situ techniques 
will require simultaneous data. 

Stratospheric NO profiles have been derived from solar absorption 
spectra recorded at sunset near 44°N [Ackerman et al., 1973; Ackerman 
et al., 1975; Louisnard et al., 1983]. These studies did not report the 
inclusion of diurnal corrections in the analysis. Our profile agrees very 
well with the results obtained in 1973 [Ackerman et al., 1973], but our 
values near 30 km are lower by * 30 percent than the measurements determined 
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Fig. 4 Comparison for NO between the present sunset infrared absorption 

measurements and in situ measurements obtained for the same latitude 
and season under high Sun conditions. 




22 


from the other flight data [Ackerman et al., 1975; Louisnard et al., 1983]. 
Measurements obtained in September 1978 with a balloon-borne pressure 
modulator radiometer [Roscoe et al., 1981] show an average daytime profile 
which decreases more rapidly than both our sunset profile and the in situ 
data plotted in figure 4. 

It needs to be emphasized that because the NO lines are weak features, 
high resolution is required for accurate quantitative measurements with the 
infrared absorption technique. Because the NO profile decreases rapidly 
below 30 km, inversion of the lower stratospheric portion of the profile is 
complicated not only by diurnal changes but also by the fact that very 
little of the NO absorption along the path originates from within the 
tangent layer. For example, for the path corresponding to the lowest Sun 
(94.14°) scan analyzed in this study, our profile indicates 77 percent of 
the NO molecules occur in layers above the tangent layer, so that small 
errors in fitting the NO lines will lead to large errors in the calculated 
tangent layer NO density. For these reasons, we estimate that the 
uncertainty in the NO profile increases from *20 percent at 32 km to 
*50 percent at 18 km. The uncertainties in the diurnal correction factors 
are not included in this estimate. 

The many factors involved in the photochemical model calculations make 
it difficult to evaluate the accuracy of the diurnal correction factors. 
Recently, we compared ground-based measurements of the total vertical column 
amount of NO with values from our photochemical calculations [Rinsland 
et al., 1984c]. The measured column amounts were derived from an analysis 
of 0.01 cm -1 resolution infrared solar absorption spectra recorded near 
sunrise and sunset on February 23, 1981, at The National Solar Observatory 
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on Kitt Peak near Tucson, Arizona. The agreement obtained in this study 
between the measured and calculated diurnal changes is close to the 
estimated precision of the measurements (=6 percent). However, it is 
important to note that the solar zenith angles of the scans studied were 
less than 86° so that the rapid changes expected for NO at sunrise and 
sunset were not observed. Measurements at a fixed altitude during sunrise 
and sunset, such as those reported by Ridley et al . [1977] and Ridley and 
Schiff [1981], provide the best test of photochemical model predictions. 

Table V presents the CO 2 volume mixing ratio values obtained from the 
analysis of the two spectral regions. The effective altitude (H e ff) 
listed for each layer is the altitude corresponding to the airmass-weighted 
effective pressure along the path through the tangent layer. The value 
listed for the layer above the balloon float altitude is the altitude 
corresponding to the effective pressure of the high sun (85.85°) scan. 

The absolute uncertainties in the C0 2 mixing ratios are estimated to be 
» 15 percent based on the uncertainties in the experimental and 
spectroscopic parameters. At this level of accuracy, our results show no 
obvious trend in C0 2 mixing ratio with altitude. Within this uncertainty, 
they also agree with the more precise (±0.5 ppmv) gas-chromatographic 
measurements made in southern France (44°N) in 1979 [Volz et al ., 1981]. 

This comparison suggests that systematic errors in the inversions arising 
from uncertainties in geometric factors (e.g., balloon float altitude, 
pointing angle) are likely to be relatively small for this data set. 

However, we note that a relatively large value of 386 ppmv was retrieved for 
the 22.3-27.4 km altitude region. 
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4.2 NO? Retrieval 

The NO 2 profile has been derived from analysis of the 1602.0 - 
1607.0 cm -1 interval (581 data points). Atmospheric lines of the v 3 and 
the v 2 + v 3 - v 2 bands of N0 2 » the v 2 bands of H 2 0, and the quadrupole 
fundamental vibration-rotation band of 0 2 are seen in this range in addition 
to several relatively weak solar CO lines [Goldman et al., 1983]. The line 
structure is superimposed on broad absorption from the collision-induced 
fundamental band of 0 2 , which extends from “ 1400-1750 cm -1 [Rinsland 
et al., 1982b]. 

Improved positions and intensities for lines of the V 3 and the v 2 + v 3 
- v 2 bands of N0 2 , derived recently from analysis of laboratory spectra 
recorded with the Kitt Peak Fourier transform interferometer, have been as- 
sumed in our analysis (R. A. Toth, JPL, private communication, 1983). For 
the N0 2 air-broadened half widths, an average value for N 2 broadening of 
0.066 cm _ 1 atm _1 at 296 K [Malathy Devi et al., 1982a] and a j -0 * 968 
temperature dependence [Malathy Devi et al., 1982b] have been assumed. Line 
parameters for H 2 0 and 0 2 were taken from the 1982 AFGL major gas line para- 
meters compliation [Rothman et al., 1983a], except for the positions of the 
stronger H 2 0 lines which were updated with the precise values measured by 
Guelachvili [1983]. No lines of CH 4 or N 2 0 have been identified within this 
interval in the stratospheric spectra [Goldman et al ., 1983]. 

The same three cases used for the NO analysis were considered. 

Table VI presents the retrieved densities, and Table VII reports the 
correction factors. Essentially all of the N0 X above 40 km is in the form 
of NO [Solomon et al., 1982], so that in computing these factors it was 
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assumed that the NO 2 contribution above the upper altitude limit of our 
photochemical model calculations (55 km) is negligible. As expected, the 
diurnal corrections are much less important for deriving the profile of N0 2 
than NO. The maximum difference in the values between the three cases is 
15 percent for the lowest layer. Since R(i,j) > 1, diurnal corrections 
decrease the retrieved N0 2 densities. Despite the improvements in the 
accuracy of the line parameters and in the analysis procedure, the profile 
obtained in the present study differs only slightly from that retrieved from 
the same set of spectra by Niple et al . [1980]. An example of the fitting 
results is presented in figure 5, where the measured and least -squares best 
fit spectra are presented for the 93.36° scan. The standard deviations of 
the residuals for the six scans range from 1.1 percent to 2.0 percent of the 
peak intensity. 

In figure 6 the results are compared with profiles derived at 33°N from 
long path infrared and visible absorption measurements at sunset by Murcray 
et al. [1974], Goldman et al . [1978], Borghi et al . [1983], and Kendall and 
Buijs [1983]. The present results are slightly lower than these published 
values. Average daytime mixing ratios derived from measurements with a 
pressure modulator radiometer near 32°N [Roscoe et al., 1981] are about a 
factor of 2 lower. We estimate that the uncertainty in our N0 2 profile 
increases from « 20 percent at 32 km to * 40 percent at 18 km. 

5. Summary and Conclusions 

We have reported simultaneous stratospheric sunset profiles of NO and 
N0 2 derived from solar absorption spectra using a nonlinear least squares 
spectral fitting technique and normalized diurnal correction factors 
calculated with a detailed time-dependent photochemical model. The NO 
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WAVENUMBER (cm -1 ) 


Fig. 5. Comparison of the 1602. 0-1 607. 0-cm‘ region of a solar absorption 
spectrum obtained at an astronomical zenith angle of 92.41° 
(tangent altitude = 27.4 km) and a least-squares best fit to the 
data. The results are presented in the same format as Fig. 2 
and Fig. 3. 



ALTITUDE (km) 




N0 2 MIXING RATIO (ppbv) 


Comparison of remote visible and infrared measurements of the NCk 
sunset profile at 31°-33°N. 
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densities obtained from the analysis are about a factor of 2 lower than 
published values obtained in situ under high Sun conditions during the same 
season and at the same latitude. This difference is consistent with the 
expected rapid decay of NO into N0 2 at sunset. The N0 2 measurements are 
slightly lower than values determined previously from solar absorption 
spectra recorded at sunset near 33°N. As a check for systematic errors in 
the analysis, the C0 2 profile has also been retrieved. These results are in 
good agreement with precise in situ measurements of C0 2 in the 
mid-stratosphere [Volz et al ., 1981]. The NO and N0 2 profiles reported here 
in combination with the simultaneous profiles of HN0 3 , H 2 0, CH4, 0 3 , and N 2 0 
reported previously [Goldman et al., 1980; Rinsland et al ., 1982a; Rinsland 
et al., 1984a] represent a consistent set of measurements useful for 
comparisons with profiles calculated by models of stratospheric chemistry. 

Our results are in agreement with previous studies [Kerr et al ., 1977; 
Murcray et al ., 1978; Boughner et al ., 1980] which show that diurnal 
corrections for NO are important for solar occultation measurements only for 
altitudes below about 30 km and that diurnal corrections produce only a 
relatively small change in the results for NO above 30 km and for N0 2 at all 
altitudes. Normalized factors calculated with the procedure described in 
this report should also prove useful for improving the accuracy of the 
retrievals of the profiles of diurnal ly varying gases from satellite 
measurements, such as will be obtained by ATMOS and HALOE. Our calculations 
indicate that normalized diurnal correction factors are relatively 
insensitive to changes in a number of model input parameters so that 
reasonably precise results can be derived without detailed knowledge of the 
profiles of the chemically-connected species, temperature, and pressure. 
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However, much further work needs to be done to access the absolute 
accuracies of the diurnal variations of NO and NO2 calculated by current 
photochemical models. Additional diurnal measurements of these species with 
simultaneous O3 and temperature data are needed, particularly with 
sufficient time resolution to study the rapid changes in concentration at 


sunrise and sunset. 
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Table I. Profile Results for NO 


Layer 

Boundaries 

(km) 

NO Sunset 

Density (10 8 cm" 

3 ) 

No Corrections 
(R=N=1) 

Case A 
Corrections 

Case B 
Corrections 

33.0-100.0 

(3.4) 

(3.4) 

(3.4) 

31.4- 33.0 

14.8 

13.2 

13.1 

27.4- 31.4 

9.1 

9.6 

9.5 

22.3- 27.4 

4.5 

6.0 

6.2 

19.7- 22.3 

3.8 

5.2 

5.7 

17.3- 19.7 

2.7 

4.5 

4.7 


Values for the layer above the balloon (33.0-100 km) are enclosed in 
parentheses to denote their larger uncertainty. 
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Table II. Correction Factors for the NO Inversion 


Layer 



R(i.j) 


N(i,j) 

Boundaries 

Scan ( ° ^ 










(km) 


Case A 

Case B 

Case A 

Case B 

33.0-100.0 

85.85 

1.0 

1.0 

1.0 

1.0 


91.30 

1.0 

1.0 

1.21 

1.21 


92.41 

1.0 

1.0 

1.12 

1.12 


93.36 

1.0 

1.0 

1.05 

1.05 


93.78 

1.0 

1.0 

1.03 

1.03 


94.14 

1.0 

1.0 

1.01 

1.01 

31.4- 33.0 

91.30 

0.98 

0.98 

1.0 

1.0 


92.41 

0.87 

0.90 

1.03 

0.995 


93.36 

0.71 

0.70 

1.03 

0.996 


93.78 

0.66 

0.66 

1.03 

0.996 


94.14 

0.63 

0.64 

1.03 

0.996 

27.4- 31.4 

92.41 

0.95 

0.95 

1.0 

1.0 


93.36 

0.74 

0.76 

1.09 

1.02 


93.78 

0.67 

0.67 

1.09 

1.02 


94.14 

0.64 

0.64 

1.10 

1.02 

22.3- 27.4 

93.36 

0.89 

0.89 

1.0 

1.0 


93.78 

0.78 

0.79 

1.12 

1.04 


94.14 

0.72 

0.73 

1.14 

1.04 

19.7- 22.3 

93.78 

0.92 

0.91 

1.0 

1.0 


94.14 

0.84 

0.82 

1.03 

1.03 

17.3- 19.7 

94.14 

0.93 

0.93 

1.0 

1.0 
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Table III. Comparison for NO between the Case B Diurnal 

Correction Factors and the Values Estimated by 
Blatherwick et al . [1980]. 


Diurnal Correction Factor R(i,j) 


Layer 

Boundaries 

(km) 

Scan 

(°) 

Case B 
(This Study) 

B1 atherwick 
et al . [1980] 

33.0-100.0 

85.85 

1.0 

1.0 


91.30 

1.0 

1.0 


92.41 

1.0 

1.0 


93.36 

1.0 

1.0 


93.78 

1.0 

1.0 


94.14 

1.0 

1.0 

31.4- 33.0 

91.30 

0.98 

1.0 


92.41 

0.90 

0.75 


93.36 

0.70 

0.50 


93.78 

0.66 

0.50 


94.14 

0.64 

0.50 

27.4- 31.4 

92.41 

0.95 

1.0 


93.36 

0.76 

0.50 


93.78 

0.67 

0.50 


94.14 

0.64 

0.50 

22.3- 27.4 

93.36 

0.89 

1.0 


93.78 

0.79 

0.50 


94.14 

0.73 

0.50 

19.7- 22.3 

93.78 

0.91 

1.0 


94.14 

0.82 

0.50 

17.3- 19.7 

94.14 

0.93 

1.0 
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Table IV. Comparison for NO Between the Case B Sunset 

Profile Derived in this Study and the Profile 
Determined by Blatherwick et al . [19803 


Layer 

Boundaries 

(km) 

NO Sunset Volume 

Mixing Ratio (ppbv) 

Case B 
(This Study) 

B1 atherwick 
et al. [1980] 

33.0-100.0 

9.36 

6.35 

31.4- 33.0 

4.40 

3.04 

27.4- 31.4 

1.95 

1.28 

22.3- 27.4 

0.59 

0.32 

19.7- 22.3 

0.31 

0.19 

17.3- 19.7 

0.17 

0.098 


The values of Blatherwick et al . [1980] correspond to the results obtained 
with photochemical effects included in the analysis. 
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Table V. C0 2 Volume Mixing Ratio 6 Obtained from 
Analysis of Two Spectral Intervals 


Layer 

(km) 

H eff 

(km) 

C0 2 Mixing Ratios 

for Spectral Range 

Average 

Value 

e(ppmv) 

1903.0-1907.0 cm- 1 
B(ppmv) 

1911.95-1917.0 cm -1 
P(ppmv) 

33.0-100.0 

37.5 

318 

299 

308 

31.4- 33.0 

31.9 

350 

345 

347 

27.4- 31.4 

28.5 

304 

341 

323 

22.3- 27.4 

23.5 

387 

385 

386 

19.7- 22.3 

20.4 

306 

328 

317 

17.3- 19.7 

18.0 

317 

340 

328 
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Table VI. Profile Results for N0 2 


Layer 

Boundaries 

(km) 

N0 2 Sunset 

Density (10 9 

cm -3 ) 

No Corrections 
(R=N=1) 

Case A 
Corrections 

Case B 
Corrections 

33.0-100.0 

(0.17) 

(0.17) 

(0.17) 

31.4- 33.0 

2.3 

2.1 

2.0 

27.4- 31-4 

2.2 

2.1 

2.1 

22.3- 27.4 

2.0 

1.8 

1.7 

19.7- 22.3 

1.1 

1.0 

1.1 

17.3- 19.7 

0.81 

0.68 

0.68 


Values for the layer above the balloon (33.0-100.0 km) are enclosed 
in parentheses to denote their larger uncertainty. 
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Table VII. Correction Factors for the NO 2 Inversion 


Layer 

Boundaries 

(km) 

Scan (°) 

R(i ,j) 

N(i,j) 

Case A 

Case B 

Case A 

Case B 

33.0-100.0 

85.85 

1.0 

1.0 

1.0 

1.0 


91.30 

1.0 

1.0 

1.58 

1.58 


92.41 

1.0 

1.0 

1.26 

1.26 


93.36 

1.0 

1.0 

1.10 

1.10 


93.78 

1.0 

1.0 

1.06 

1.06 


94.14 

1.0 

1.0 

1.02 

1.02 

31.4- 33.0 

91.30 

1.01 

1.02 

1.0 

1.0 


92.41 

1.06 

1.06 

1.00 

0.971 


93.36 

1.14 

1.18 

1.00 

0.969 


93.78 

1.16 

1.20 

1.00 

0.969 


94.14 

1.17 

1.21 

1.01 

0.969 

27.4- 31.4 

92.41 

1.02 

1.03 

1.0 

1.0 


93.36 

1.10 

1.12 

1.01 

0.998 


93.78 

1.12 

1.16 

1.02 

0.999 


94.14 

1.12 

1.17 

1.02 

1.00 

22.3- 27.4 

93.36 

1.03 

1.10 

1.0 

1.0 


93.78 

1.06 

1.09 

1.10 

1.02 


94.14 

1.07 

1.12 

1.11 

1.02 

19.7- 22.3 

93.78 

1.02 

1.04 

1.0 

1.0 


94.14 

1.04 

1.08 

1.05 

1.06 

17.3- 19.7 

94.14 

1.03 

1.06 

1.0 

1.0 




2. Government Accession No. 


3. Recipient's Catalog No. 


4 . Title and subtiti^t ra t 0S ph e ri C NO and NOo Profiles at Sunset 
from Analysis of High-Resolution BaTloon-borne Infrared 
Solar Absorption Spectra obtained at 33°N and Calculations 
with a Time-Dependent Photochemical Model 


7. Author(s) 

C. P. Rinsland, R. 
A. Goldman**, F. J. 


E. Boughner, J. C. Larsen*, 
Murcray**, and D. G. Murcray** 


9. Performing Organization Name and Address 

NASA Langley Research Center 
Hampton, VA 23665 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D. C. 20546 


15. Supplementary Notes 

* Systems and Applied Sciences Corporation 
Hampton, VA 23665 


5. Report Date 

August 1984 


6. Performing Organization Code 

147-44-02-70 


8. Performing Organization Report No. 


10. Work Unit No. 


11. Contract or Grant No. 


13. Type of Report and Period Covered 

Technical Memorandum 


14. Sponsoring Agency Code 


** University of Denver 
Denver, CO 80208 


16. Abstract 


Simultaneous stratospheric vertical profiles of NO and NO 2 at sunset have 
been derived from an analysis of 0.02-cm -1 resolution infrared solar 
absorption spectra recorded from a float altitude of 33 km with the 
University of Denver interferometer system during a balloon flight from 
Alamogordo, N.M., near 33°N on October 10, 1979. A nonlinear least squares 
procedure has been used to analyze the spectral data in regions of 
absorption by NO and-N0£ lines. Normalized factors, determined from 
calculations of time-dependent altitude profiles with a detailed 
photochemical model, have been included in the onion-peeling analysis to 
correct for the rapid diurnal changes in NO and NO 2 concentrations with time 
near sunset. The C0 2 profile has also been derived from the analysis and is 
reported . 


17. Key Words (Suggested by Author(s)) 

Remote sensing 
Atmospheric optics 
Infrared spectra 
Atmospheric gases 
Nitric oxide 


19. Security Classif. (of this report) 

Unclassified 


18. Distribution Statement 


Unclassified — Unlimited 
Subject catagory 46 


20. Security Classif. (of this page) 

21. No. of Pages 

Unclassified 

45 


22. Price 

A0 3 


N-305 


For sale by the National Technical Information Service, Springfield, Virginia 221S1 



















